Probing Majorana fermions in spin-orbit coupled atomic Fermi gases 
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We examine theoretically the visualization of Majorana fermions in a two-dimensional trapped 
ultracold atomic Fermi gas with spin-orbit coupling. By increasing an external Zeeman field, the 
trapped gas transits from non-topological to topological superfluid, via a mixed phase in which both 
types of superfluids coexist. We show that the zero-energy Majorana fermion, supported by the 
topological superfluid and localized at the vortex core, may be visible through (i) the core density 
and (ii) the local density of states, which are readily measurable in experiment. We present a 
realistic estimate on experimental parameters for ultracold *'^K atoms. 
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Majorana fermion [l| - particle that is its own antipar- 
ticle - has attracted considerable attentions from a wide 
area of physics A particular interest comes from its 
non-Abelian exchange statistics that is crucial for topo- 
logical quantum computation 0| . Two well separated 
Majorana fermions may form a non-local fermionic state, 
as a non-local qubit for inherently fault-tolerant quantum 
memory. As a portal for future quantum technology, the 
realization of Majorana fermions in a highly controllable 
manner is of great importance and a timely quest. 

Majorana fermions are believed to exist in a number 
of two-dimensional (2D) physical settings, including frac- 
tional quantum Hall states at filling i/ = 5/2 \a, vortex 
states of px + ipy superconductors/superfluids [1,0|, and 
surfaces of 3D topological insulators in proximity to a 
s-wave superconductor Q. Majorana fermions may also 
emerge in ID quantum systems with strong spin-orbit 
coupling, such as quantum wires Q and optically trapped 
ID fermionic atoms [l^. AH these appealing proposals 
are yet to be realized experimentally. 

In this work, we examine the possibility of observing 
Majorana fermions in the vortex core of a spin-orbit cou- 
pled ultracold atomic Fermi gas in 2D harmonic traps. 
This is a scenario discussed earlier by several researchers 
based on a theoretical concept originated from 
Jackiw and Rossi, who predicted the existence of vortex- 
core Majorana fermions in (2+1) dimensional Dirac the- 
ory [l^ . Here, we perform a fully microscopic calcula- 
tion with Bogoliubov-de Gennes (BdG) equation, which 
enables simulations with realistic experimental parame- 
ters. Our study is motivated by the recent creation of 
non-Abelian gauge fields in a Bose-Einstein condensate 
(BEG) of *^Rb atom s llSl and its possible realization in 
fermionic "'"K atoms |16| . 

We find that by increasing a Zeeman field, a topological 
superfluid emerges from the trap edge and extends grad- 
ually to the whole Fermi cloud, supporting zero-energy 
states (ZES) at the vortex core and trap edge. This topo- 
logical phase transition is detectable through a sudden 
change of the atomic density inside the vortex core, as- 
sociated with the occupation of the Majorana state. We 



show that the wave function of the Majorana fermions 
can be inferred from the local density of states (LDOS) 
at the core. 

Mean-field BdG equation. — We consider a trapped 
2D atomic Fermi gas subject to Rashba spin-orbit cou- 
pling Vso(r) = —iX{dy+idx) and a Zeeman field h, which 
may be prepared in a single pancake-like optical trap 
V{r,z) = Mluj'^r'^ + uj'^z'^]/2 with trapping frequencies 
ujz ^ uj±. We note that 2D Fermi gas has recently been 
realized in experiments (TtI. [l8|. The system is described 
by H = fdr [Ho {r)+Hi (r)] , where 

Ho(r) = ''^i^^(i-)^- + [p\Vso{r)i^i + H.c] (1) 

and Hi{r) ~ C/o'0j(r)V-'J(r)V't(r)V't('') describes the con- 
tact interaction between opposite spins. Here ^ are the 
creation field operators for the spin-up and -down atoms, 
= -h'^V^/{2M) + Muj\r'^/2 ~ fi-ha^ is the single- 
particle Hamiltonian in reference to the chemical poten- 
tial fi. The interaction strength Uq is to be regularized 
via l/C/o + Ekl/(ft'kVA/-H^a) = M/{47:H^)lniEa/E). 
Here Ea is the binding energy of the two-body bound 
state (lol . [20| and i? > is the relative collision energy. 

The low-energy fermionic quasiparticles are solved by 
the mean-field BdG approach, 'HBciG^'r; (r) = -E'r)^'?; (r). 
Using the convention for Nambu spinors '5^ (r) = 
[u-fri,u^ff,v-f^,v^jj]'^ , the BdG Hamiltonian reads. 
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is the order parameter, to be solved self-consistently 
in conjunction with the atomic densities, rig- (r) = 

(l/2)E,[k-.|V(i?r,) + \v^^\' fi-E^)]. Here f (x) ^ 
1/ (e^/''^-^ -f l) is the Fermi distribution function. 
The chemical potential ^ is determined by the to- 
tal atom number N = J dr[n^ (r) 4- n^(r)]. With 
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Figure 1: (color online) (a) Phase diagram, along with the lowest eigenenergy of Bogoliubov spectrum, (b) Configuration of the 
mixed phase ai h — 0.5Ef- The local critical field is calculated in the vortex-free state, (c) The gap profile with and without 
vortex at /i = 0.5Ef- (d), (e), and (f) Energy spectrum at h/Ep = 0.4, 0.5, and 0.6 (where Ef = tfk%/{2M) = VNhuj± 
is the Fermi energy) in the presence of a single vortex. The color of symbols indicates the mean radius \/{r^ /rp (where 
rp = {4:Ny/'^^yh/{Mu}±) is the Fermi radius) of the eigenstate, which is defined by (r^) = fr^[\u^f + \uif + \v^f + \vif]dr. 
The color of symbols changes from blue when the excited state is localized at the trap center to red when its mean radius 
approaches the Thomas- Fermi radius. The CdGM states are indicated by blue squares in (d) and (e), and by a solid line in 
(f). The dashed lines show the inner edge and outer edge branches. 



a single vortex at trap center, we take A(r) = 
A(r)e~"'' and decouple the BdG equation into differ- 
ent angular momentum channels indexed by an inte- 
ger m. The quasiparticle wave functions take the form, 
[ut^(r)e-^'^,?x^^(r),«t^(r)e»'^,^;^^(r)]e'■(™+l)'^/^/2^. We 
have solved self-consistently the BdG equations using 
the basis expansion method. For the results presented 
here, we have taken N = 400 and T = 0. We have used 
Ea = 0.2Ep and Xkp/Ep = 1, which are typical param- 
eters that can be readily realized in a 2D ""^"K Fermi gas 
[isj . Other sets of parameters, with varying interaction 
strength, SO coupling strength and temperature, have 
also been tried. 

The use of Nambu spinor representation leads to an 
inherent redundancy built into the BdG Hamiltonian. 
'H-BdG is invariant under the particle-hole transformation, 
Wcr (r) — )■ V* (r) and iJ^ — >■ —E^. Thus, every eigenstate 
with energy E has a partner at —E. These two states 
describe the same physical degrees of freedom, as the 
Bogoliubov quasiparticle operators associated with them 
satisfy F^; = tLb- Iii the expressions for order parameter 
and atomic density, this redundancy has been removed 
by multiplying a factor of 1/2. 

Majorana fermions. — The particle-hole redundancy, 
however, is very useful to illustrate a non-trivial fea- 
ture when the Zeeman field h is beyond a threshold, 
h > he = a//?+A2 and the system is in a topologi- 
cal state Q, hosting ZES within the energy gap. Due 
to i? = 0, the associated quasiparticle operators satisfy 
Tq — Fq. Thus, a zero-energy quasiparticle is its own 



antiparticle - exactly the defining feature of a Majorana 
fermion Because of the redundant particle-hole rep- 
resentation, the ZES or Majorana fermion is a /la// of or- 
dinary fermion and thus must always come in pairs. Each 
of the paired states, localized separately in real space, can 
hardly be pushed away from i? — by a local perturba- 
tion [j, 0| , giving rise to the intrinsic topological stability 
enjoyed by Majorana fermions. It is straightforward to 
show from the BdG Hamiltonian that the wave function 
of Majorana fermions should satisfy either Ua (r) = v* (r) 
or u„ (r) = —V* (r). The former follows directly from the 
particle-hole symmetry. For the latter, the related quasi- 
particle operator satisfies Tq = ^ f J. This is necessary 
for expressing an orindary fermion using paired Majorana 
fermions (2l|. 

Phase diagram. — Figure [T] reports the phase dia- 
gram (a) along with the quasiparticle energy spectrum 
of different phases (d, e, and f) in the presence of a sin- 
gle vortex. By increasing the Zeeman field, the system 
evolves from a non-topological state (NS) to a topologi- 
cal state (TS), through an intermediate mixed phase in 
which NS and TS coexist. The mixed phase, unique for a 
trapped system, can be easily understood from the point 
of view of local density approximation, in which the lo- 
cal chemical potential /i(r) = n — Ma;^r^/2 and the or- 
der parameter A(r) decrease continuously away from the 
trap center. As a result, the local critical Zeeman field 
hc{r) = •\/^2(r) -I- A2(r) becomes smaller at the trap 
edge, as shown in Fig. (lb) for h = Q.5Ep. This cre- 
ates a ring of TS at the outer region where h > hc{r), 
surrounding the inner region which is non-topological. 
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The full topological transition occurs when the Zeeman 
field is larger than the critical field at the trap center, 
h > -^//i^ + A2(0), where A(0) is the gap at trap center 
in the absence of the vortex. 

The topological phase transition into TS is well char- 
acterized by the low-lying quasiparticle spectrum, which 
has the particle-hole symmetry Em+i = — -E_(m-i-i)- As 
shown in Fig. (fd), the spectrum of the NS is gapped 
by order parameter, even at the trap edge. In the mixed 
phase (Fig. (le)), however, three branches with small 
energy spacing appear. The branches labeled as "Outer 
edge" and "Inner edge" consist of the eigenstates with 
wave functions localized at the edges of the TS as An- 
dreev bound states, as indicated by the arrows in Fig. 
(lb). The two edge states at to = —1 have nearly zero 
energy. The other branch, shown by blue squares, is a 
series of discrete localized states at the vortex core, i.e., 
the so-called Caroli-de Gennes-Matricon (CdGM) states 
[2^ . When the TS extends over the whole cloud (Fig. 
(If)), the dispersion of the inner edge branch moves into 
the continuum. The energy of the state with to = — 1 in 
both the CdGM branch and outer edge branch becomes 
essentially zero. 
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Figure 2: (color online) Wave functions of the two lowest- 
energy modes in the full topological state aX h = Q.&Ef, (a) 
and (b), and in the partial topological state at ft = 0.5Ef, (c) 
and (d). The wave functions are in units of a/ A4uj± / h . 

The eigenstates with nearly zero energy at to = —I, 
i.e., the two edge states in the mixed phase as well as 
the outer edge state and CdGM state in the full TS, can 
be identified as zero-energy Majoranan fermions in the 
thermodynamic limit. To show this, we plot in Fig.[2]the 
wave function of these eigenstates in the full TS (a and 
b) and in the mixed phase (c and d). In both cases, we 
observe a bond and anti-bond hybridization between two 
well-localized wave functions: one satisfies (r) = v* (r) 
and the other Ua (r) = — w* (r), which is exactly the 
symmetry of the wave function required by Majorana 
fermions. The hybridization is caused by a quasiparti- 
cle tunneling between paired Majorana states [23, which 



leads to the splitting of degenerate zero-energy states to 
finite energies ±i?zES- The tunneling barrier between 
the two edge states is lower, so the energy splitting is 
relatively larger (i.e., i?zES ~ 0.006£^f)- In contrast, the 
tunneling between the outer edge state and CdGM state 
seems to be more difhcult, giving rise to an exponentially 
small splitting (i.e., i?zES ~ 4.8 x 10~^Ef). 

The different phases can therefore be identified from 
the lowest eigenenergy of energy spectrum, as plotted in 
Fig. (la). Within the NS, it decreases slowly as the Zee- 
man field increases. The decrease becomes exponentially 
fast with the appearance of a partial TS and a minimum 
is reached when the cloud just becomes fully TS. Further 
increase of the Zeeman field will reduce the order pa- 
rameter and hence the tunneling barrier between the two 
Majorana states, which are localized respectively at the 
vortex core and trap edge, leading to a steady increase of 
the lowest eigenenergy. We note that the exponentially 
small energy of Majorana fermions in the full TS, inher- 
ent to the finiteness of trapped Fermi cloud, should be 
suppressed by increasing the number of total atoms. 

Probing Majorana fermions. — In the TS, the occu- 
pation of the Majorana vortex-core state affects signifi- 
cantly the atomic density and LDOS of the Fermi cloud 
near the trap center, which in turn gives an unambiguous 
experimental signature for observing Majorana fermions. 




Figure 3: (color online) Zeeman field dependence of spin- 
up (a) and spin-down (b) densities at the vortex core. The 
density is normalized by the Thomas- Fermi density np = 
{■\/N /tt) ^Mojx/ft- The insets show the core density distribu- 
tions at ft = O.G-Ef- The red dot-dashed lines show the result 
by excluding artificially the Majorana vortex core state, whose 
contribution is shown by the shaded area. 

Figure [3] presents the spin-up and -down densities at 
the trap center, n^ifi) and n^(0), as a function of the 
Zeeman field. In general, n^(0) and 714^(0) increases and 
decreases respectively with increasing field. However, we 
find a sharp increase of n^(jS) when the system evolves 
from the mixed phase to the full TS. Accordingly, a 
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change of slope or kink appears in ri-t-(O). The increase of 
n^(0) is associated with the gradual formation of the Ma- 
jorana vortex-core mode, whose occupation contributes 
notably to atomic density due to the large amplitude of 
its localized wave function. We plot in the inset of Fig. 
(3b) ni{0) at /i = 0.6i?_F, with or without the contribu- 
tion of the Majorana mode, which is highlighted by the 
shaded area. This contribution is apparently absent in 
the NS. Thus, a sharp increase of n^{0), detectable in 
in-situ absorption imaging, signals the topological phase 
transition and the appearance of the Majorana vortex- 
core mode. This feature persists at typical experimental 
temperature, i.e., T = O.lTp. We note that experimen- 
tally it is more favorable to take a time-of-flight imaging 
of the cloud after an expansion time, in order to have 
enough resolution to visualize the vortex core. We antic- 
ipate the sharp increase in ^^^(O) may persist for a short 
expansion time. Alternatively, we may tune quickly an 
external magnetic field to take the Fermi system to the 
BEC limit across Feshbach resonances. In this way, the 
vortex core can be imaged clearly after the time-of-flight 
just like in an atomic BEC. 
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Figure 4: (color online) Log-scale contour plot of the LDOS 
for spin- up (a) and spin-down (b) atoms at ft = Q.QEp- Here 
we use r = 0.01-Bf. The LDOS is units of uf/Ef. 

In the full TS, the wave function of the Majorana 
mode can be determined by measuring LDOS through 



spatially resolved rf-spectroscopy (23 . which pro- 
vides a cold-atom analog of the widely used scanning 
tunneling microscope in solid state. We show in Fig. |4] 
the spin-up and -down LDOS at ft, = O.QEp defined as 
p,{v,E) = {l/2)Y.^[\u„,\' 5{E-E^) + \v„,,\' 5{E + E^)l 
where the (5-function can be simulated by a Lorentzian 
distribution with a suitable energy broadening F. In- 
side the vortex core, the contribution from the Majorana 
mode and other CdGM states is clearly visible within 
the superfluid gap. In the case of F, fc^T < AS, where 
/^E - A2(0)/(2-Bf) is the energy spacing of CdGM 
states [12 which for typical parameters as used in our 
calculation turns out to be about 10 nK, the Majorana 
fermion contribution Pcr,zEs(r,0) may be singled out. As 
P(T,ZEs(r, 0) cx |w^^(r)|^ = |t;^^(r)|^, the spatially re- 
solved rf-spectroscopy maps out directly the wave func- 
tion of the Majorana vortex-core state. 

In closing, we note that 2D ultracold atomic Fermi 
gases are an ideal platform for probing and manipulating 
Majorana fermions because of the unprecedented con- 
trollability and flexibility. This is particularly useful for 
the purpose of topological quantum computation, using 
Majorana fermions as qubits For instance, two 2D 
atomic Fermi gases formed by a double-well potential 
along z-axis, each of which has a single vortex at the cen- 
ter, can host four Majorana fermions for carrying out the 
basic information process. In this configuration, inter- 
well quantum tunneling of Majorana fermions is possible, 
providing another potential means to detect the interest- 
ing topological phase transition. 
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